Decoherence Control in Open Quantum System via Classical Feedback 



O 
O 

> 
O 

(N 



Narayan GanesarQ and Tzyh-Jong TarrQ 

Electrical and Systems Engineering. 
Washington University in St. Louis 

In this work we propose a novel strategy using techniques from systems theory to completely elim- 
inate decoherence and also provide conditions under which it can be done so. A novel construction 
employing an auxiliary system, the bait, which is instrumental to decoupling the system from the 
environment is presented. Our approach to decoherence control in contrast to other approaches in 
the literature involves the bilinear input affine model of quantum control system which lends itself 
to various techniques from classical control theory, but with non-trivial modifications to the quan- 
tum regime. The elegance of this approach yields interesting results on open loop decouplability 
and Decoherence Free Subspaces(DFS). Additionally, the feedback control of decoherence may be 
related to disturbance decoupling for classical input affine systems, which entails careful application 
of the methods by avoiding all the quantum mechanical pitfalls. In the process of calculating a 
suitable feedback the system has to be restructured due to its tensorial nature of interaction with 
the environment, which is unique to quantum systems. The results obtained are qualitatively dif- 
ferent and superior to the ones obtained via master equations. Finally, a methodology to synthesize 
feedback parameters itself is given, that technology permitting, could be implemented for practical 
2-qubit systems to perform decoherence free Quantum Computing. 
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I. INTRODUCTION 



Various authors have studied control of decoherence 
of an open quantum system. Decoherence Free Sub- 
spaces(DFS) help preserve quantum information in an 
open quantum system. However, the presence of symme- 
try breaking perturbations or control hamiltonians acting 
on an open quantum system which is essential to perform- 
ing arbitrary transforms in the system hilbert space 7i s , 
could also lead to loss of information by inevitable trans- 
fer of states out of DFS, due to the nature of the control 
hamiltonians. Hence this renders the quantum system 
at best a noiseless memory, much less a dynamic quan- 
tum computer, whose state needs to be transformed in 
order to perform computations. Recently Lidar and Wu 
Kielpinski et. al.jH, Brown et. al have 
proposed a combination of open loop bang-bang pulses, 
universal control in order to perform computation within 
the DFS via control pulses. In this work we propose a 
novel strategy, exploiting the geometry of the bilinear 
control system on the analytic manifold to completely 
eliminate decoherence in the presence of symmetry break- 
ing control hamiltonians and still preserve complete con- 
trollability of the system in order to perform arbitrary 
transforms. We also explore the possibilities and provide 
conditions under which it can be done so. This unified 
approach to control of decoherence lets us analyze the 
open loop decoupling problem which directly leads us to 
the existence of DFS and secondly closed loop decoupling 
via a classical feedback to the control system which leads 



us to robust decoherence control. This work is a continu- 
ation of the previous results [3] wherein some of the theo- 
retical groundwork was laid to study the problem of open 
loop decoupling, which are now extended to closed loop 
control and feedback design here. The approach used 
here is fundamentally different from approaches adopted 
by other authors in that (z) the bilinear form of control 
system is used which is amenable to classical systems the- 
oretical results instead of the stochastic master equation 
for the state evolution, (ii) the approach does not aim at 
mitigating or slowing down the decoherence rate rather 
aims at completely eliminating via a suitable non-linear 
feedback. The experimental feasibility is discussed for 
a finite state environment acting on a two qubit system 
which is a rather reasonable approximation. A procedure 
to compute the feedback using the invariant subspace for 
a system is provided. A detailed step by step algorithm 
to determine the invariant subspace itself on the tangent 
space T^(M.) is also provided. In order to compute the 
feedback parameters a good estimate of state of the sys- 
tem is essential. A reliable information extraction scheme 
utilizing indirect continuous measurement via a quantum 
probe in the context of a decohering quantum system was 
studied inllal. 



II. PREVIOUS WORK 

Consider an open quantum system interacting with the 
environment described by, 

= [H ® T e (t, x) +l s ® H e (t, x) + H SE (t, x) 
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Here the argument x denotes the spatial dependance of 
the combined system-environment state x) as well as 
control hamiltonians Hi, and where Ui are the strength 
of the control respectively. Ho, He, H$e are the system, 
environment and interaction hamiltonian acting on TL Sl 
Ti e and TLs®7~i e (system, environment and the combined) 
Hilbert spaces respectively. For ease of notation we will 
suppress the spatial dependance. Define an output equa- 
tion which could either be a non-demolition measurement 
or a general bilinear form given by, 



Theorem II. 1. Let 

Co = C(t) 

C n = span{ad J Hi C n -i(t)\j = 0, 1, . . . ; i = 1, . . . , r} 
C n =\ (ads + ^Y C n ;j = 0,1,- •■ 1 



y(t) = m\C(tM(t)) 



(1) 



where again C(t, x) is assumed to be time- varying oper- 
ator acting on the system Hilbert space. For instance for 
a finite system the non-hermitian operator C — \m){n\ 
when plugged in eq. ([1]) would yield the coherence be- 
tween the respective states of the system or for an electro- 
optic system the operator C — aexp(iujt) + (V exp(-iujt) 
would yield the output of a real non-demolition observa- 
tion performed on the system. In order to study the in- 
variance properties with respect to the system dynamics 
of the above time dependent quantum system, we define 
f(t,x,ui, ■ ■ ■ ,u r ,H S B) = y(t,0 for t £ [to,*/] to be a 
complex scalar map as a function of the control functions 
and the interaction Hamiltonian Hsb over a prescribed 
time interval. The function / is said to be invariant or 
the signal y(t, £) is said to decoupled from the interaction 
Hamiltonian Hsb if, 



f(t,X,Ux, • 



-,H S b) = f(t, x,u x , ■ ■ ■ ,u r ,0) (2) 



for all admissible control functions m, • • • ,u r and a given 
interaction Hamiltonian Hsb- Then the condition for 
such an output signal to be decoupled from the interac- 
tion hamiltonian in the open loop case is given by the fol- 
lowing theorem [LJ|, which follows an iterative construc- 
tion in terms of system operators. 



The vector fields Ko = 



Ki = 





Hi£(x, t) 



K p = 



1 


H S pt{x,t) 



and K] = 







corresponding to drift, control and inter- 



H SB i{x,t) 

action can be identified to contribute to the dynamical 
evolution. It was already noted that the the system was 
said to be decoupled if it satisfied equations [31 namely, 



Define a distribution of quantum operators, C(t) = 
A{C 1 {t),C 2 (t), ■ ■ -,C n (t), ■■■}. The output equation (??) 
of the quantum system is decoupled from the environ- 
mental interactions if and only if, 
Case (I): Open Loop, 



[C(t),HsB(t)} = 



(4) 



Case (II): Whereas the necessary conditions for Closed 
Loop control is, 

[C,H SB }=0 
[C(t),H SB (t)](ZC(t) 

In this work we will be primarily concerned with de- 
signing feedback for quantum systems of the form u = 
a(£) + P(S,)v where a and f3 are real vector and a full 
rank real matrix of the state (or its estimate thereof) of 
dimension 1 x r and r xr respectively. We examine a few 
systems of interest with control hamiltonians, that might 
be decoupled via feedback of the above form. 

Definition II. 1. The vector field K T satisfying equa- 
tions (0i is said to be in the orthogonal subspace of the 
observation space spanned by the one-forms 

dy(t, £),dL Kio y(t, £),•• • , dL KiQ ■■■ L Kin y(t, £),•• • (5) 



V0 < i 0> • 



< r and n > 



Denoted by K T £ O x 

Lemma II. 2. The distribution O x is invariant with re- 
spect to the vector fields Ko, ■ ■ ■ , K r under the Lie bracket 
operation, (i.e) if K T G O , then [K T ,Ki] £ O for 
i = 0, • • • , r 



III. A SINGLE QUBIT SYSTEM 



Recalling, 



Lk,L 



Kj^k zo ■ ■ ■ 



L Kl y(t,£) = 

L Ktn y( t >O = 



(3) 



Consider a single qubit spin-1/2 system coupled to a 
bath of infinite harmonic oscillators through an inter- 
action hamiltonian Hsb- The hamiltonian of the sys- 
tem+bath can be written as, 



H 



LOO 



k 



u k b\b k 



■9kh) 
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FIG. 1: The 2 Qubit system is allowed to interact with an- 
other qubit, the bait whose interaction with the thermal bath 
is controlled. 



where the system is acted upon by the free hamiltonian 
Hq and the decoherence hamiltonian Hsb- As is well 
known there is a rapid destruction of coherence between 
|0) and |1) according to the decoherence function given 
by [HI • In order to cast the above problem in the present 
framework we consider a bilinear form of an operator C 
that monitors coherence between the basis states. Con- 
sidering C to be the non-hermitian operator |0)(1| we 
have a function y(t) given by y(t) = (£(t)|C|£(i)) that 
monitors coherence between the states |0) and The 
problem now reduces to analyzing the applicability of the 
theorem III. 1 1 to the given system. It can be seen right 
away that the condition [C, Hsb] for the distribution 
C defined previously, as calculated in the previously 
This implies that the coherence is not preserved under 
free dynamics or in presence of open loop control. In 
order to eliminate this decoherence by feedback we now 
assume the system to be acted upon by suitable con- 
trol hamiltonians {Hi, • • • , H r } and corresponding con- 
trol functions {ui, ■ ■ ■ , u r }. As we pointed out earlier 
the necessary condition is relaxed to [C, Hsb] C C, with 
the operators C and Hsb still required to commute with 
each other [C,Hsb] = 0. For the single qubit example 
the second condition fails to hold, again as outlined [lj], 
thus leaving the system unable to be completely decou- 
pled and hence vulnerable to decoherence even in the 
presence of closed loop and feedback control. 



IV. TWO QUBIT CASE 

In case of two or multiple qubits there always exist De- 
coherence Free Subspaces (DFS) that are immune to the 
decohering hamiltonian. Recently, Fortunato et. al[29| . 
Mohseni et. al[3(|, Ollerenshaw et. al[3l| proposed and 
demonstrated computation within the DFS. However it is 
not certain that the system could be contained within the 
DFS at all times under the action of the control hamil- 
tonians a x , a y for the system. With the effort to steer 
within the DFS the authors of above work could show an 
improvement to previous methods, but still prone to ef- 
fects of decoherence. As a simple calculation suggests 



that with the initial state ci|01) + C2 1 10) , within the 
DFS for a 2-qubit system and after a time t of control 
acting on the first qubit which transforms the state to 
ci(cos<|0) +sint|l))|l) +c 2 (cosi|l) +sint|0))|0) which is 
clearly out of the DFS. Recently Lidar and Wu [11], [13, 
Kielpinski et. al.[13], Brown et. al [33| have proposed 
a combination of open loop bang-bang pulses, univer- 
sal control and DFS in the context of ion trap quan- 
tum computers to perform computation within the DFS 
via control pulses, which again produces an improvement 
over previous results but still prone to decohering effects. 
However we follow a different control strategy where in 
we seek to completely eliminate the influence of Hsb 
based on feedback control and a novel construction in 
order to perform decoherence free control. The corre- 
sponding 2-qubit control system can be written as, 



dim 

dt 



2 z 



5> fc &£&J |£(t)) (6) 



+ ( Ul (t)*W +U2 (t)o-y + U3 (t)*W + u 4 (t)a^Mt)) 

which satisfies the basic necessary condition 
[C, Hsb] = but not the stronger condition pro- 
vided in Case(w) of the theorem. Hence the system 
would eventually leave the DFS and is susceptible to 
decoherence in the presence of arbitrary control, in other 
words, not entirely decoupled from Hsb- In order to 
analyze the system and the conditions in the presence 
of a classical state feedback u = a(£(t)) + (3(£(t)).v 
the corresponding conditions (it) of the theorem are to 
be examined. Since the operator Hsb € B(H a ® Tie), 
the set of skew hermitian linear operators acting 
non-trivially on both system and environment hilbert 
space, whereas the operators in the distribution C for 
the above control system is confined to B(TL S ) that act 
trivially on the environment hilbert space. Hence the 
necessary condition specified in Theorer dII.il would not 
be satisfied non-trivially unless the distribution C acted 
non-trivially on both 7i 8 and 7i e . In other words the 
distribution includes operators of the form ^ A a ® B a 
for a countable index set {a} and operators A a and B a 
operating on system and environment respectively. The 
above forms cannot be achieved by control hamiltonians 
acting only on the system. However the situation can 
be salvaged if one considered a "bait" qubit whose 
rate of decoherence or the environmental interaction 
can be modulated externally at will and the bait qubit 
is now allowed to interact with our qubits of interest 
through an Ising type coupling. With the help of the 
following construction we will be able to generate vector 
fields of the form Ki artificially, which will be seen 
to provide great advantage. With the coherence func- 
tional y(t) = (£(i)|01)(10|£(i)) where the state 
vector is now the total wave function of system+ bait+ 
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environment. Both the qubit systems are assumed to 
interact with the same environment with the additional 
requirement that the bait qubit's decoherence rate be 
controllable. Physically this amounts to a coherent 
qubit with controllable environmental interaction. The 



scalability and advantages of this construction are 
analyzed in the next section. 

The Schrodinger equation for the above system can 
now be written as, 



J 



!B 3f> 



+u 3 (t)a x ^ + u 4 (t)aW + ?£ a W + u 5(7 W + + U7 J ia ^a^ + u 8 J 2 a^a^ + u 9 ]T a z h \w k b\ + w* k b k ) £(*)> 



r 



where cr x ,a y ,a z are regular hermitian operators and 
u\ (t) to ug (t) are time-dependent piecewise constant con- 
trol functions. The terms of controls 117 and u 8 are gener- 
ated by the Ising type coupling between qubits 1, 2 and 
the bait with the corresponding coupling constants J\ 
and J 2 respectively. The last term in the above control 
system is due to the interaction of the bait qubit with 
the environment whose interaction enters the system in 
a controllable way hence can be treated as a separate 
control hamiltonian. Keeping in mind the following com- 
mutation relations between different pairs of operators, 



2ia z 



>y,< 



[h,b\,] = S k k 



= 2ia x 
[bk,blb k ] = b k 



2icr v 



[blblb k ] 



through control hamiltonians alone. But in order to gen- 
erate such a vector field one has to know before hand 
and as time progresses the exact values of the environ- 
mental coupling coefficients g k which at best could only 
be described by a stochastic process. Hence in the light 
of the aforementioned difficulty just rendering the co- 
herence independent of Hsb seems like a much better 
alternative. 



a z = |l)(l|- 10X01,^ = |0)<1| + |1><0| 

<7 v =i|0)<l|-i|l><0| 

and C = |01)(10| = (a x 1} - ia ( y ] ) ® (ct£ 2) + ^ 2) )/4, 
we have [C, Hsb] — and [C,Hsb] for instance con- 
tains terms of the form a x ® ®^{g k b\ + g k b k ) which 
are not zero. Fortunately with the above construction 
these terms can be seen to be present in the distribution 
C, which can obtained under the sequence of operations 
[C,H{\ = cwzQcTy, [[C,Hi],H 5 ] = c 2 <J z ®<J x ®J2(.9kbi + 
g* k b k ), [[[C, Hi\,H h \,H 2 ] = c 3 a x ® 1^ ® J2(g k bl + g* k b k ) 
and the corresponding <j y term is obtained via the se- 
quence, [[[C, H 2 ], H5], Hi]. Since both terms are present 
in C, so is their linear combination. Hence both the nec- 
essary conditions as outlined by the theorem for closed 
loop decouplability are satisfied for the above system. 
Hence we are one step closer to decoupling the coher- 
ence between the qubits from Hsb- In fact it can be 
seen that the operator Hsb itself can be generated by 
the control hamiltonians through the lie bracket opera- 
tion H SB = [[H 5 ,H 2 ],Hi] or [[H 5 , Hi], H 2 ]._ Hence any 
term in [C,Hsb] is trivially contained in C. Hence, it 
might seem at first that the effects of Hsb on the sys- 
tem could be nullified by generating an equivalent —Hsb 



V. SCALABILITY 



It can also be seen that the above approach works for 
finite number of qubits coupled to only one bait qubit 
through the same a x ^a y ^ interactions. Such an interac- 
tion can be implemented using the same technology nec- 
essary for multi-qubit quantum computers wherein a fi- 
nite number of qubits are entangled to a single qubit that 
is capable of readout and storage of an oracle's query re- 
sults. With the underlying theory of disturbance decou- 
pling in place all that remains now is synthesis of the feed- 
back control itself. Since the conditions [C, Hsb] C C 
and [C,Hsb] = turn out to be necessary conditions, 
with the proof of sufficiency requiring further insight into 
design and construction of appropriate control fields we 
will for the next few sections follow an alternative for- 
malism called an Invariant Subspace which is a part of 
the tangent space T^(Ai) of the analytic manifold. It will 
be seen later that the two seemingly different approaches 
viz. (i) the conditions in terms of operators of the system 
and (ii) The tangent space formalism, complement one 
another in terms of obtaining a complete solution to the 
problem of disturbance decoupling. 
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FIG. 2: The isobar of y(f) is represented by the sphere and 
the nullspace ker(dy(t) is a tangent to the sphere at the point 

VI. INVARIANT SUBSPACE FORMALISM 

Consider the necessary and sufficient conditions for de- 
couplability 

L Kl y(t) = (8) 
L Kl L Ko y(t) = (9) 

Hence LkqLk^^) = 0. The above equations af- 
ter subtraction imply L^ KoK .^y(t) — 0. The other 
necessary conditions viz. £[Ko jK^KiVit) = and 
LKjL^K^t) = imply that L [[K(uK ^ K]] y(t) = 0. In 
fact the above pattern of equations could be extended to 
any number of finite lie brackets to conclude that 

L[[-[K ,K Z1 IK I2 ]-K lk ]y(t) = (10) 

which leads us to define a set of vector fields or distribu- 
tion A that share the same property, 

K v G A s.t L Kv y{t) = (11) 

It is observed immediately that Kj E A. Such a distri- 
bution A is said to belong to null space of the function 
y(S,,t). And from the necessary conditions listed above 
the distribution is observed to be invariant under the con- 
trol and drift vector fields Kq, ■ ■ ■ , K m , (i.e) \IK V G A, 

[K u ,Ki\ G A,Vz G 0,- • • ,m 

Simply stated, 

[A, Ki] c A, Vi G 0, • • • , m (12) 

We will henceforth refer to the distribution as the in- 
variant distribution. It is also to be noted that the above 
calculations are reversible and the original necessary and 
sufficient conditions can be derived starting from the in- 
variant distribution. Hence the necessary and sufficient 
conditions for open loop decouplability can now be re- 
stated in terms of the invariant distribution. 



Theorem VI. 1. The output y(t) is unaffected by the 
interaction vector field Ki if and only if there exists a 
distribution A with the following properties, 

(i) A is invariant under the vector fields Kq, Ki, ■ ■ ■ , K m 

(ii) Ki G A C ker(dy(t)) 

Hence existence of the invariant subspace is essential 
to decouplability of the system in question. It is now all 
the more important to determine the invariant subspace 
(if any) for the given system and output equation. In 
order to compute the invariant distribution it properties 
discussed above comes in handy and provides a means to 
go about computing the distribution as well. 

The procedure starts out by assigning the entire null 
space ker(dy(t)) to invariant distribution A and succes- 
sively removing parts of the distribution that don't sat- 
isfy the other properties (i.e), invariance with respect the 
vector fields Kq, ■ ■ ■ , K r . In other words, remove parts of 
A whose lie brackets with Kq, ■ ■ ■ , K r do not lie within 
A. Of course, the above mentioned procedure involves 
computing inverse image of Lie brackets as described be- 
low. 



A. Invariant Distribution Algorithm 

Algorithm 1: 
Step 1: Let A = ker(dy(t, £)). 

Step 2: A i+ i = A, - {6 £ A t : [S, Kj] $ A 4 < j < r} 
Step 3: Maximal invariant distribution is such that 
A* = A, when A l = A i+i , Vi . 

The above is an iterative procedure that computes dis- 
tributions Ai in order to arrive at the final invariant dis- 
tribution A* = A. Where the '— ' is the set removal 
operation. Let us redefine the set to be removed as, 

Si = {5e A t : [5,Kj] £ A 4 ,V0 < j < r} 

Hence the set Si can also be written as, 

Si = invQA,-,!^] - A,), V0 < j < r (13) 

where 'inv' is the set theoretical inverse mapping of the 
linear map [., Kj],VQ < j < r, taking values in A i.e, 

inv(r) = {5 G A; : [S, Kj] = t, V0 < j < r} 

Figure([3]) outlines the schematic of the algorithm. 

One of the foremost issues to be addressed is the con- 
vergence of the algorithm. However at this point we are 
not fully equipped to study the converge as the proof 
below will introduce additional ideas to discuss conver- 
gence. It is to be noted here that A^ is always a distri- 
bution^ vector space) for all i. Hence the set Si is such 
that, the removal of Si from A^ results in a distribution 
of lower dimension Aj_|_i. Hence removal of the set Si 
removes a subspace A^ contained within the distribution 
Ai. Hence we have A i+i c Aj and 

A i+1 + Ai = A, 
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{6 : [6, Ki) £ A} {8: [6, K % \ e A} 



FIG. 3: Shaded portions (light and dark) mark the original 
distribution A* C ker(dj/(t)). The dark shaded portion repre- 
sents the core of the distribution that is invariant and the light 
shaded portion, the part of distribution that is not invariant 
and the white portion, image of [., Ki] that lies outside Aj. 

Where the '+' now denotes the direct sum of two sub- 
spaces. However, The procedure outlined above is not 
convenient as it involves solving for inverse mapping un- 
der lie bracket operation [131 It is for this reason that we 
would like to perform the calculations in the orthogonal 
complement within the dual space T£ (M) of the tangent 
space. The algorithm can now be reformulated entirely 
in terms of the orthogonal complement, ilj C T£(M) of 
the distribution Aj. (i.e) the inner product, 

(w, 5) = 0, S Q and 5 £ A 

denoted by (fi, A) = or f2 = A- 1 . Hence the algo- 
rithm now starts out by setting f^o = span(<iy(i)) and 
iteratively adding the subspace that was removed by the 
previous removal operation and finally inverting the co- 
distribution so obtained to recover A, (i.e) 
Step 1: Set fi = span(c?y(<, £)). 
Step 2: = Qi+ (Aj)*. 

Step 3: The Algorithm converges to 12* = fit when 
fij+i = flj,Vi 

where (.)*(not to be confused with A*) stands for the 
corresponding dual vectors within the dual space T£(A4), 

(i.e), if Aj = span{5i, ■ • • 5k} then 

(Aj) = spanjc^! , • • • where (u/j,#j) = 1. Now the 

task at is to determine the subspace (Aj) . It is helpful 
to examine the relations between the distributions Aj, Aj 
and Aj+i. Note that, dim(Aj) + dim(f2j) = N and Aj+i 
is orthogonal to Aj. In fact it can also to be seen that 
Aj+i is precisely, 

A i+1 = {S£ A, : [5, Kj] £ Aj, V0 < j < r} (14) 



in order to locate the subspace Aj we have to determine 
the complementary subspace (look for vectors that are or- 
thogonal) to eq. (fT4|) and within Aj. From the identities 
of Lie derivatives, 

L K] (u, 6) = {L Kj u, S) + (u, [6, K t ]) (15) 

Hence for uj £ f2, and 6 £ Aj+i C Aj, we have, (u>, S) — 
and {w,[S,Kj\) = 0(eq.([I31)). Hence (L Kj lu,S) = 
from the previous identity (|15|) . In other words Lx^i is 

orthogonal to Aj+i. Since A i+ i is orthogonal to (A)j 
and f2j, we have, 

L Kj Qi COj + (Aj)* (16) 

Now consider the same equation (fl~5|) . for all S £ Aj and 
all uj £ f2j we have (w, 5) = 0. But since [S, Kj] ^ Aj we 
have (oj, [5, Kj]) ^ for some ui £ Hi, hence (Lk^, 6) ^ 
as well. Hence for any S £ Aj there exists an u> € Oj 
such that (I^Wjif) ^ 0. (i.e) 

(Aj)* C L Kj SU (17) 
Hence from eq. (fTo| and (jTTJ) we conclude that 

(Aj)* + fi< = Qj + Lif 3 .f2j,VQ <j<r. (18) 

although it is possible to prove the stronger condition, 
f2j + A* = LKjQi- We now state the algorithm without 
proof: 

Step 1: Set fig — sp&n(dy(t, £)). 
Step 2: = Qi + Ljr (Qj) + ££ =1 W^O- 
Step 3: The Algorithm converges to fl* — fij when 
= Qj, Vi. 

Maximal invariant distribution A is such that A* = 
rr- 1 . As seen in the proof each step of Algorithm 1, re- 
moves a set from a vector which amounts to removing 
a finite dimension limited by dimension of the tangent 
space. Hence the convergence of the algorithm is depen- 
dent on the finite dimensionality of the tangent space at 
point £(t) which can be guaranteed by the finiteness of 
control Lie Algebra, which will be studied in the following 
sections. 



B. Observation space and Tangent Space 

In definition III.ll the observation space spanned by 
dy(t,0,dL Kio y(t,0,---, dL KiQ ■ ■ ■ L Kin y(t,£), •••VO < 
*0> • • " j in < t and n > was defined and it can be eas- 
ily seen that the necessary and sufficient condition for 
open loop decouplability ([3]) is equivalent to being or- 
thogonal to the observation space according to def. III.ll 
The orthogonality relation also follows from the simple 
Lie derivative identity, 



which is a restatement of Step 2: of Algorithm 1 . Hence 



L Kr y(t,0 = (K T ,dy(t,0) (19) 
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From[14j, it can be seen that the one forms 
dy(t, £),dL Ki y(t, £) etc can be expressed in terms of the 
commutators of operators and hamiltonians, G, i?o,-ffi- 
Infact the operations performed in the observation space 
provide an alternative formulation to the theory devel- 
oped in terms of the tangent space and invariant dis- 
tributions. As can be seen the structure of the out- 
put equation y(t) = (£{t)\C (t)\£(t)) made possible the 
simplifications of Lie derivatives of scalar functions to 
commutators of operators and enjoys ease of calculations 
when compared to computing Lie derivatives of vector 
and co- vector fields, if one were to compute the invari- 
ant subspace. Hence it is to be noted that the necessary 
and sufficient conditions for open loop decouplability can 
just be stated in terms of the observation space without 
ever having to calculate the invariant distribution which 
is precisely what Theorem III. II sets out to do. And it is 
also to be noted that the Theorem is a consequence of 
the orthogonality relation in the observation space (Def- 
inition III . 1 1) . 

However when it comes to feedback decouplability the 
two different formalisms play equally important roles 
in constructing a quantum system that might be de- 
coupled using feedback. The observation space formal- 
ism provides important necessary conditions (in terms of 
the commutators of operators) while designing a quan- 
tum control system while the tangent space formalism 
is indispensable to calculating the feedback parameters 
a(£(t)), f3(£(t)) once the system of interest is known to 
be decouplable using feedback. 

VII. SYNTHESIS OF FEEDBACK 
PARAMETERS a(£),/3(£) 

In this section we study the explicit formulation of the 
feedback control that ensures complete decoupling of the 
coherence functional from Hsb- It is to be seen that 
this formulation can be applied to outputs other than 
the coherence functional we wish to monitor, like that of 
a non-demolition observable. 

Definition VII. 1. A distribution is said to controlled 
invariant on the analytic manifold D u if there exists a 
feedback pair (a, [3), a, vector valued and (3, matrix val- 
ued functions such that 

[£ O ,A](OCA(0 (20) 
[Xi,A](0cA(0 (21) 

where, 

r 

K = Ko + ^ajKj 

3=1 

and 

r 



It is to be noted that Kq and Ki are the new drift and 
control vector fields of the control system after applica- 
tion of feedback (a, (3). The problem of decoupling via 
feedback can now be cast in the original framework of 
open loop decouplability by requiring that the feedback 
vector fields now satisfy the open loop decouplability con- 
ditions viz. 

[£o,A](Oc A(0 
[#i,A](0c A(£) 

and that A be contained entirely within the null space of 
the output function (i.e), 

A c kei(dy) 

With the above characterization of feedback decoupla- 
bility the task now reduces to finding a distribution that 
might satisfy the above invariance conditions with re- 
spect to the feedback vector fields, (Kq, K%, ■ ■ ■ ,K r ), 
which in turn requires the knowledge of the feedback 
functions a and (3. What seems to be a deadlock sit- 
uation can now be resolved by further simplifying the 
invariance condition stated above. 

Lemma VII. 1. An involutive distribution A defined on 
the analytic manifold T> u is invariant with respect to the 
closed loop vector fields (Ko, K\, • • • , K r ) for some suit- 
able feedback parameters a(£) and /3(£) if and only if, 

[X ,A]cA + G (22) 
[K it A] c A + G (23) 

Where G is the distribution created by the control vec- 
tor fields. 

G= span {#!,••• ,K r } (24) 

At this point it is possible to express the necessary 
and sufficient conditions for the feedback control system 
(K , Ki, ■ ■ ■ ,K r ) to be decoupled from the interaction 
vector field Kj just as we were able to provide conditions 
for open loop decouplability. Moreover the conditions 
can be expressed entirely in terms of the open loop vec- 
tor fields and the controlled invariant distribution with- 
out ever having to involve the feedback parameters a(£) 
and The following theorem provides the conditions, 

Theorem VII.2. The output y(t,g) = (£|C(t)|£) can be 
decoupled from interaction vector field Kj via suitable 
feedback (a, [3) if and only if there exists an involutive 
distribution A such that, 

[Kq, A] C A + G 
[Ki,A] c A + G 

and A C ker(dy) 
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Proof. ( ) The following proof covers the lemma as 
well as the theorem above. Assuming that A is locally 
controlled invariant or in other words invariant with re- 
spect to the closed loop vector fields (K , K\, ■ ■ ■ , K r ) for 
some feedback parameters and within an open 
set in V^j. If r £ A, then it can be seen that, 

r r 

[K u r] = \^K h r\ = £ [Kj , r] - £(L T /3y)*i 

3=1 3=1 

as we know the left hand side is still contained within A 
and the last term on the right side is a linear combination 
of vectors that generate G. Hence 



It will be seen that the above requirement will be easily 
satisfied for the extension to control algebra to be dis- 
cussed following this proof. It is also to be noted that 
we haven't imposed any non-singularity restrictions on 
the distributions above. Now the control vector fields 
Ki € G could be written as a linear combination of the 
vector fields {|t>i), • • • , \vn)} at each point £. 



d N 

Ki = J2 c v \ V J) + c v \ v i)> VI < » < r 

3=1 3=d+l 

K = Kf + K°{ where K? G A and K° <£ A}. 



J2^[K,,r}eA + G 

3=1 

and since (3 is assumed to be nonsingular it is possible 
to solve for individual [Kj , r] by mere inversion of the 
matrix fcj and can be found to be linear combination of 
vectors in A + G and hence, 

[K u t] e A + G 

Now consider, 

[K ,r] = [K + 1 £ aj K j ,T] 

r r 

= [K , r]+Y^ a 3 {K 3 ,T] - ^{Lra^Kj 

3=1 3=1 

Since the left hand side belongs to A and since 
[Kj,r] e A + G, 1 < j < r it can be immediately seen 
that [K Q , t] e A + G as well. 

(<=-) For the proof of sufficiency the following geo- 
metric visualization is helpful. Let the dimension of the 
distribution A be d. Since A is involutive there exist 
d vectors fields, locally non-vanishing in a neighborhood 
C/C^n^of f , {| Vl ), • • • , \v d )} E T S (M) that are 
linearly independent and, 



A = span{|«i), • • • , \v d }} 



(25) 



s.t [\vi), \vj)] € A, VI < i,j < d. Now let the dimen- 
sion of A + G at £ be d + q. It is now possible to 
find another q linearly independent vector fields labeled 
{\v d+1 ), • • • , \v d+q )}, such that [\vi), \vj)] £ A,V1 < i < 
d, d + 1 < j < d + q. As a special case one could think 
of a local co-ordinate basis that are mutually commut- 
ing and linearly independent. Let the dimension of the 
tangent space at the point £ be N. Finally it is possible 
to find N — d — q additional linearly independent vectors 
that complete the vector space T^(M), by Gram-Schmidt 
procedure or otherwise (i.e), 



The vector fields are devoid of components in A. And 
since dimension of A + G is d + q it can be seen that 
the r vectors K°, • • • , K° span a q dimensional subspace. 
Hence it is always possible to generate q linearly inde- 
pendent vectors and r — q zero vectors via suitable linear 
combinations of K°, • • • , K°. Let the linear combinations 
be such that, 



r N 

Y,(3i 3 k; = \v d+1 )+ ~°io\ v o) 

3=1 j=d+q+l 
r N 

Y,NK° 3 = \v d+2 )+ Yl 52 >j> 

3=1 j=d+q+l 



r N 

j2p qj K° = \ Vd+q )+ 5 «m 

3=1 ]=d+q+l 

and 



T^(M) = span{>i), • • • , \v d ) \v d+1 ), 
\vd+ q +i), ■■■ , \vn)} 



\Vd+q), 



(26) 



3=1 



3=1 



The flij matrix so formed is precisely the feedback pa- 
rameter that is used to generate the closed loop vector 
fields Ki, 1 < i < r. 



Ki = PijKj denoted by (3.K 



In order to prove this we note the action of the 
above linear combination on the open loop vector fields 
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dim(A) = d 



FIG. 4: The dimension of controlled invariant distribution is 
d and the control distribution G is partitioned into {Kf} and 
{K° }. The basis vectors \vi), ■ ■ ■ , \vd) span A. 



Ki, ■ ■ ■ ,K r ,(i.e), 



r d N 

Y $ x i K i = Y ^ + + Y 5l J 1^ 

r d JV 

Y K 3 = Y I V 3 ) + I V d+2 ) + Y ^ I V 3 ) 
3=1 3 = 1 j=d+q+l 



and for 1 < i < q, 

d N 

[r, (3 l3 Kj } = [ T ,Y cij \vj) + \v d+i ) + Y ~ c i3 1 v 3 ) 

3=1 3=d+q+l 
d N 



By noting that r 6 A, A is involutive and \vi), • • • , \vd) 
commute with \vd+i), ■ ■ ■ , |ujv), [\vi)> \ v j}] € A,V1 < i < 
d, d+l<j<d + q the above equation can be seen to 
simplify to, 

r N 

[T,Y^3Kj}eA+ Y 

3 = 1 j=d+q+l 

but since we already have [r, /Sylfj] G A + G from 1(25]) 
the above relation is possible only if i T (cy ) = 0. Hence 
we have, 

r 

J=l 

The argument is trivial for q + 1 < i < r and it can be 
easily seen that [r, 0ijKj] G A for all 1 < i < r and 
r G A. Now in order to construct the feedback parameter 
a, by an argument analogous to (|29[) . we can show that, 

r 

[r, # ] = [r, #o + Y a i R i\ e A + G 



r d N 

Y Pl3 K 3 = Y C lj\ V 3) + \ y d+q) + Y Cc u \ V 3 ) (27) 
3 = 1 3 = 1 j=d+q+l 

and 

r d 

Y Ah-i,^ = Y + 

3=1 3 = 1 



Y^ K ^Y^\ v ^ + 

3=1 3=1 

where the first terms on the right hand side of the above 
equations can be seen to be from (3K d and the later terms 
from (3K° . We will suppress the summation for ease of 
notation and all the following terms below are assumed 
to be summations from 1, • • • , r in the recurring index 
variable. Now from the necessary conditions we have, 

[t, Kj] G A + G, Vt G A and 1 < j < r. (28) 

and hence, 

[r, PijKj] = 0ij [r, Kj] + L T (fi^K, G A + G. (29) 



because both [r, K ] and [r, Ki] belong to A + G. Let, 

d d+q N 

K a = Y^ J \v J )+ Y c o\ v j)+ Y c >i> ( 3 °) 

3=1 3=d+l j=d+q+l 

It is now possible to find a suitable linear combination of 
right hand side of equation set (|2T|) and the above ([50)) 
in order to form Kq, 

d d+q N 

K = Ko+otiKi = Y^\ v i)+ Y k 3\ v i)+ Y 5 i\ v i} 

3=1 3=d+l j=d+q+l 

where fcj's are constants w.r.t £ and t where as c,-'s are 
some functions of £(t). In particular by a suitable linear 
combination, fcj's can all be made zero. It can again be 
seen that for all r G A, 

N 

[r,K ]eA+ Y M^)M 

j=d+q+l 

and hence L T (cj) are equal to zero in order to satisfy the 
necessary conditions and hence, 

[t,K ] G A 

The closed loop drift vector field was formed by setting 
Kq equal to Kq + a./3.K for a suitable row vector a. 
Hence the feedback parameter a = a. (3. □ 
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In addition to proving the necessary and sufficient 
conditions we have also outlined a procedure to compute 
the feedback parameters ck(£) and /?(£) from the maximal 
controllability invariant distribution A, which elicits 
the application of Tangent space formalism in output 
decoupling. Hence it is imperative that we compute the 
maximal invariant distribution A for the synthesis of 
feedback. From the necessary and sufficient conditions 
we see that the distribution A has to satisfy conditions 
(|20I21|1 or equivalently (|22H23[) and that A C ker(dy) 
for complete decouplability. Obviously (|22M23p has 
the advantage that we do not need the knowledge of 
feedback parameters. Now, similar to the open loop 
case we can formulate an algorithm in order to arrive at 
the much sought after invariant distribution, the general 
idea being: Start out by assigning the whole of null 
space of y(t) to A and iteratively remove the part of the 
distribution that does not satisfy conditions (|22H23p . 



Step 1: Let A = ker(dy(i, £))■ 

Step 2: A l+1 = A, - {i5 e A, : [6,Kj] £ A, + G,V0 < 
3 < r} 

Step 3: Maximal invariant distribution is such that 
A* = Ai+i = A,. 

Employing the same logic as before in determining 
the open loop invariant distribution we can perform the 
computation in the dual space TZ(M) and arrive at the 
following algorithm which is easier to compute, 

Step 1: Let Q = span(dt/(t, £)). 

Step 2: Q i+1 = n i + L JfD (fi i nG- L )+ELiijf 4 (^nG J -). 

Step 3: The Algorithm converges to VL* = fii+i = 



VIII. EXTENSION TO CONTROL ALGEBRA 

In the previous sections we provided a state feedback 
given by the vector a(£) and matrix which were as- 
sumed to be analytical functions of the state £. In partic- 
ular, the analyticity is required for the proof of necessity 
as well as sufficient conditions. However, the class of ana- 
lytic functions is too restrictive in terms of feedback that 
can actually be implemented on the system. For exam- 
ple, by rapid pulses which are arbitrarily strong and fast 
one can generate lie bracket of the vector control vector 
fields which can act as a new control to the system avail- 
able for feedback. In the light of non-analytic feedback it 
might be necessary to modify the conditions that guar- 
antee decouplability of the system. Another approach 
which is sufficiently general would be to use the theory 
already developed for analytic feedback to systems whose 
control vector fields belong to the control algebra of the 
original system, (i.e) we propose to use the system, where 
Ki G {Ki, • • • , K t } L a = Q- The theory of analytic feed- 
back can now be extended to controls from the control 
algebra instead of just the original set of controls. Hence 



we can restate the conditions for decouplability in terms 
of the control algebra, which follows directly from the 
previous theorem as, 

Lemma VIII. 1. The output y(t) is decouplable via an- 
alytic feedback functions a(£) and /?(£) from the interac- 
tion vector field Ki if and only if there exists a control- 
lability invariant distribution A, (i.e) 



[A,g]cA©6 
[A,C] cAfflg 



(31) 
(32) 



where C — {ad J K ,Ko,i 

{#!,••■ ,K r } LA l 



1, 



,r;j = 0, 1 • • • } and Q 



The above lemma just states a condition and does not 
provide an explicit formulation of the application of feed- 
back. In order to provide the analytic feedback we con- 
sider a modified system with additional control vector 
fields generated from the original system. Consider the 
following modified system with finite dimensional control 
algebra Q, 



dt 



= K \m)+J2^i\m+KM(t)) (33) 



where the vector fields A; 6 Q which are generated by 
the vector fields of the original system are such that Q = 
span{A' l5 • • • , A m }, (i.e) the set of vector fields K i: not 
necessary a linearly independent set form a vector space 
basis for Q. This is a required condition as the analytic 
feedback functions which can only generate utmost linear 
combinations of the existing control vector fields, (i.e) 
span{Ai, • • • , K r } is inadequate to leverage the set of all 
possible controls. Hence it is necessary to modify the 
original system in order to utilize the repertoire of all 
possible controls for efficient feedback control. It is also 
to be noted that in so doing we do not alter the set of 
reachable or controllable set of the original system, but 
altering the output decouplability instead which is an 
observability property of the system. 



IX. EXAMPLES 

As an example of the above formalism consider a single 
qubit and a two qubit system coupled to the environment, 

= y<^(*) +J2 uJkb l bk ^ + u ^m + u 2 <j v m 



J2 <7 ^kb{ + g* k b k )m 



with the output, 



y(t) = (m\c\m) 



When we check against the necessary condition, 
J2k a z(9kbl + 9kh)^(t) e ker(dy(i)) which we notice the 
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single qubit system fails to satisfy, the conclusion that a 
single qubit system is not decouplable coincides with re- 
sults obtained earlier by operator algebra. Now, consider 
the following two-qubit system eq.© 



It can be clearly seen that the interaction vector field 
in deed belongs to A/ = Ej, fe ^(Sfc&I + 9k b k)£(t) 6 
ker(dy(t), where j = 0, 1 and k = 0, 1, • • • , but 



dt 



LOq 



+E (E^J (gkbt+g* k b k )\m) 

+ (ui(t)aW +u 2 (t)o-« + u 3 (t)a^ + )|£( t )) 

Which has a DFS of dimension 2, span{|01), |10)}, the 
states within which remain coherent in the absence of 
controls. The real problem arises in the presence of sym- 
metry breaking perturbations or control hamiltonians. 
Hence the problem at hand is to render the states coher- 
ent even in the presence of arbitrary control. Consider 
the output of the form, 

y (t) = m\c\m) 



j 

= c - E rfiidkbl + g* k b k )\0, eg. for ,i = 1, 2 



up to a constant c, neither belongs to the span of the 
control vector fields, control algebra generated by the 
above vector fields or the controllability invariant distri- 
bution A. The last condition can be seen by the fact that 
[Ki, Ki] does not belong to ker(dy(£)) and hence does not 
belong to A C kei(dy(t)) either. Now consider the two 
qubit system with bait, which was discussed in the ear- 
lier section. The control system governing the mechanics 
following the Schrodinger eq.© is given by, 



mi 

dt 



E Y*® + E ^ b l b * «*) + E E a ^ G& 6 1 + 9th) ffl + { Ul {t)a^ + u 2 (t)aW + u 3 (t)a^ 
vJ=i k J k \ j J 

+n 4 (i)af + ^af) + u 5 a^ + u 6 a^ + ujj^a^ + u s J 2 a^a^) + u 9 £ a^{w k b\ + w* k b k )i{t) 

k 

(34) 

I 



with o~ x \ y \ z iiow skew hermitian and the same output 
equation as before. It is seen that Kj £ ker(dy(t)) and 



[KM = [a<J>6 J2 <W (9kbl + 9th)t] 
j 

= c -^yi(9kb{ + gih)\0 



now belongs to the control algebra generated by the addi- 
tional vector fields introduced by the bait system. Hence 
the system which was designed in order to meet the nec- 
essary condition, [C,Hsb] C C, given by the observa- 
tion space formalism is also seen to meet the conditions 
given by tangent space or controllability invariant dis- 
tribution formalism. A rather interesting scenario arises 
when the drift vector field Ao is a part of the ideal of Q 
and the interaction vector field Kj which is a part of the 
invariant subspace A C ker(dy(t)), is already contained 
within the control algebra, (i.e) Kj £ Q. The necessary 
and sufficient conditions for decouplability using feed- 
back are trivially satisfied as [Kj, JQ] £ QNK £ Q and 



[A/, A ] £ Q. Hence, 



[A,Ki] C A: 

[A, Ko] C A * 



Q 
Q 



(35) 
(36) 



and the invariant subspace A can now be guaranteed to 
exist and at least one dimensional equal to span{Ao}. 
Hence existence of feedback and decouplability is guar- 
anteed for the above system. 



X. THE CONTROL SYSTEM 

In the previous section we had only discussed a brief 
outline of the implementation of disturbance decoupling 
for quantum systems. In this section we present the con- 
struction of actual control system and the control vector 
fields. The bait qubit as discussed before was primarily 
used to get a handle on the environment so we may gen- 
erate vector fields that could help decouple the system 
from the vector field Kj. Let the following denote the 
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various hamiltonians acting on the system, 

j=l k 
HSB= ^ fe^J (fffc6 * +ff * 6fc) 

H 9 = J2^ b) (^kbl + wp k ) (37) 

k 

and let us denote by K[, the vector fields generated by 
the hamiltonian ifj, (i.e), Ki = Hi\i/j). Now consider the 
particular back and forth maneuver via controls uq and 

u%{t) = 1, and ug(r) = 0, for r S [0,t] 
t/6(t) = 0, and Uq(t) = 1, for r S [t, 2t] 
U6(r) = — 1, and «g(r) = 0, for r e [2t, 3t] 
U6(r) = 0, and ug(r) = —1, for r £ [3i, 4t] 

The corresponding unitary time evolution operator at the 
end of time instant At is given by, 

UUt) = e (- lH zt) e (-iH<A) e (iH 6 t) e (iH<A) 

= exp(- i [H 6 , H 9 ]t 2 + 0(t 3 )) 

the series expansion by Campbcll-Baker-Hausdorff for- 
mula. In the limit that t = dt — > 0. The effective di- 
rection of evolution is given by the commutator of the 
corresponding hamiltonians, but to the second order in 
time. Hence we could devise a control vector field in the 
direction given by the commutators of the corresponding 
hamiltonians Hq and Hg, where, 

[He,H 9 ]=c.ai b ^(w k bl+w* k b k ) 

k 

where c is a real constant for a skew hermitian Hq and 
Hg . In fact it is possible to generate any direction of evo- 
lution with arbitrary strength corresponding to repeated 
commutators of the hamiltonians Hi ■ ■ ■ Hg of the phys- 
ical system (1341) . In order to compute commutators of 
tensor product operators we use the following identitiy, 

[A ® B, C <g> D] = CA <g> [B, D] + [A, C] <gs BD 

With another control field H% entering the picture we 
could generate the following direction in conjunction with 



the previous maneuver [[H S ,H 5 ], [H 6 ,Hg]], 

= c 1 .[J 2 *<?ty b \*®^(w k bt + w* k b k )] 

k 

= c.a^a^J2M{ + wlb k ) (38) 

k 

Consider the similar maneuver between controls 
i*4, uq and us, which generates the direction of evolution 
corresponding to the following repeated commutator, 

[H 4 ,H 8 } = [a^,J 2 a^a^} = c.a^a^ (39) 

where c is a real constant for a skew hermitian Hi,H^. 
Again, from operating on equations (|38)) and (|39[) we get, 

k 

k 

= c.aV.l( b \J2(^ b i + w l b K) (4°) 

k 

Hence we have generated an effective coupling between 
qubit 2 and the environment with the help of the bait 
qubit and its interaction with the environment and qubit 
2. It is important to note that the hamiltonian so ob- 
tained by the above control maneuver now acts trivially 
on the hilbert space of the bait qubit, a property which 
will be found to be extremely useful later. It is also possi- 

(2) 

ble to generate the cri counterpart of the above coupling 
by a similar maneuver, given by, 

c. CT i 2 >.lW.£K&£+<6 fe ) (41) 

k 

Again by a symmetric and totally similar argument we 
could generate a coupling between the environment and 
qubit 1, which would be given by, 

c.a^l^.J2( w k b l+ w *k b k) and (42) 
k 

c.aW.l^.^kbi + wlb,) (43) 
k 

Now noting that the constants c in the above equations 
could be controlled independently and arbitrarily, we can 
write the preliminary form of the actual control system 
which achieves disturbance decoupling. Gathering terms 
(j40| - (|43|) . we construct the following control system for 

mm 

dt ' 
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E + E ^ b l b * + E E (»4 + + + u ^v ] + + u *< 2) m 

\3=1 k J k y j J 

E(^ 6 I + w *m + u ^v ] E( wfe6 I + w * b *) + E^ fc6 I + w k h k) + m «4 2) J2( wkb i + w *m J £(*) 

V fc fc fe k / 

(44) 

I 



By restructuring the control vector fields as above we are 
hoping to capture the entire control algebra by a simple 
linear span of the control vector fields which is essential to 
analytical feedback theory. Let us again, investigate the 
dccouplability of the above control system from its neces- 
sary conditions, that {i)Kj £ Ac ker(dy), (ii)[Ki, Ki\ G 
A + G, where G — span(iTi • • • , K$), is the distribution 
generated by the control vector fields above. By consid- 
ering [K^Kx] = oi^EfcOfc&l + w%b k )\£{t)) , which is 
already contained within G. The conditions are also sat- 
isfied for the vector fields K^, K 3 and K 4 . However with 
the vector field K 5 , we note that [Kj, K 5 ], 



EJE^J (9kbi+g* k b k )m, 

k 



where w.l.o.g w k = C\.g k for an arbitrary constant c\ G C. 
For an infinite dimensional environment the vector fields 
that contain higher powers of ^2 k (w k b\. + w k b k ), cannot 
be expressed as a linear combination of its lower powers 
as can be seen from its action on a particular number 
state \n), 



{w k b\ + w* k b k )\n) 
{w k b\ + wlb k ) 2 \n) 



=w k y/n\n — 1) + w* k \Jn + l|n + 1) 
--2\w k \ 2 n\n) + wly/n(n - l)\n - 2) 
+ wfy/(n + l)(n + 2)\n + 2) 



for some n. In other words the above term is neither con- 
tained in G nor in A C kcr(dy), because L r y(£_,t) ^ 0. 
The only way to correct the above situation is to in- 
clude the vector r as a control vector field in the con- 
trol system above. This can be achieved by similar 
maneuvers between the vector fields above, (i.e), r — 
c.[Ki, [Ki, #2£]]- Now again, since r is a new control 
vector field, it must satisfy condition (ii) above. But 

[if/, r] = c.aP (j2 k ( w kb{ + w£&fc)) £(t), now generates 
the next higher power of the same environmental term, 
which necessitates us to find a way to include that in 



our control vector fields as well. In fact, it is possible 
to generate any power of the environmental term by re- 
peated commutators, which is linearly independent of all 
the previous terms and hence generates a new direction 
of flow within the analytic manifold. And it is impossible 
to include all the successive powers in our control vector 
fields. Hence the best we could hope to achieve under the 
present circumstance is to obtain an approximate solu- 
tion to disturbance decoupling. It is to be noted that the 
above problem arises only in an infinite dimensional en- 
vironment and restricting the dimension of environment 
is a reasonably good approximation. Hence we present a 
experimentally realizable scheme to demonstrate the the- 
ory of disturbance decoupling to practical quantum sys- 
tems. The following system captures the essence of the 
problem as well as the solution itself. Before we present 
the example we summarize the results obtained thus far 
in a concise form. The following table is helpful in noting 
the above decouplability results, 





Open 
Loop 


Closed 
Loop 


Closed Loop 
Restructured 


Single Qubit 


NO 


NO 


NO 


Two Qubit 


NO 


NO 


NO 


Two Qubit or higher 
with bait qubit 


NO 


NO 


YES* 



*-The system can be completely decoupled under the 
additional assumption of a finite dimensional environ- 
ment. 

We note that the conditions for dccouplability from 
Open loop to Closed loop to Closed Loop Restructured 
are progressively relaxed. Hence a system that is not 
Closed Loop Restructured decouplable cannot be Closed 
Loop or Open Loop decoupled. 

Finite State Environment Environment always ap- 
pears to be in a stationary state (also called the Gibbs 
State). An essential element of the stationary state which 
is most stable and extremely resilient is the coherent 
state of an electromagnetic system. Coherent states is 
generated by the action of the displacement operator 
D(a) = e ( aa1 - a " a ) on the vacuum state |0). An elec- 
tromagnetic system when perturbed from one coherent 
state simply settles in another coherent state. It is la- 
beled by a complex number a, that denotes the strength 
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of the state. The state is given by, 

00 n 

|a) = e -W£ ° I 

^"1 vn! 

n=u 



(45) 



where |n) is the number state. It can be the seen that the 
coefficients of higher n decrease rapidly and since squared 
sum of the coefficients is convergent, with major contri- 
bution from lower states it is a reasonable approximation 
to neglect higher energy states of the electromagnetic 
system. In fact, this is the basis for the experimental 
realization of "dual rail optical photon quantum gates" , 
where in only the |0) and |1) photon states are used to 
represent the system under the premise that contribu- 
tions from higher energy photons are negligible. Hence 
we consider the following model for a finite state har- 
monic oscillator with N energy states which will be later 
dubbed as the environment. The creation and annihila- 
tion operators act on the system as follows, 



a\n) = \Jn + l\n + 1) for n < N 
a\n) = for n > N 

a* \n) = \pn\n — 1) for n < N and n > 
a'\n) = for n > N and n = 



(46) 



It was recently shown by Fu et. al.[18[ in their model of 
truncated harmonic oscillator that such a system up to 
energy state iV was feasible. Hence the schematic pre- 
sented here can be readily implemented if one were able 
to create and sustain a controllable interaction between 
the electromagnetic and spin system(the bait). Now con- 
sider a single spin- 1/2 system with hamiltonians a z , a x , 
and a y . The state of the system is represented as tp = 
[cq ci]' in the vector form where the coefficients corre- 
spond to the two states and cq, c\ € C s.t |co| 2 + |ci| 2 = 1. 
The tangent vector to the system is given by the action 
of the skew hermitian operators on the state ip, (i.e), 

i> = Vx\y\z'4 ! = [ c 'o c 'i}' where c ,Ci G C. In other 
words, in order to express any vector in the tangent 
space as a real linear combination of other vectors we 
require at least 4 linearly independent vectors given by 
o'zV'? o'xip, cr y ip, hp. For the case of a 2 spin-1/2 system 
the number of linearly independent vectors required is 
8, given by a subset of <Ji ® <Jjip for {i,j} € {x,y,z,Q\. 



For the case of 2 spin-1/2 system coupled to a 3 state 
environment, the tangent space is 4 x 3 dimensional and 
the number of linearly independent vectors required to 
span the entire tangent space is 4 x 3 x 2 = 24. In other 
words we require 24 linearly independent control vector 
fields to make absolutely sure that the conditions for de- 
couplability are met. Let the environment be governed 
by a single 3 level harmonic oscillator. The different en- 
ergy levels are given by {|0), |1), |2)}, and a general state 
in this basis is given by |e) = cq|0) + ci|l) + c 2 |2). We 
can now examine the linearly independent vectors gener- 
ated by the powers of the bath/environment operator for 
the three level system by taking into account the defin- 
ing relations (|4"6"]l . As it can be seen that the following 6 
vectors, 



IV = co|0)+ci|l)+c 2 |2) 
{wb^ + w*b)ip = 

|0) + W2c 2 |l) + w*Co\l) + w *V2 Cl \2) 
(wtf + w*b) 2 %p = 

[w 2 (b ] ) 2 + ww*tfb + w*wbtf + w* 2 b 2 ]ijj 

(wtf + w*b) 3 ip = [w 3 ^) 2 + w 2 w*(b ] ) 2 b + w 2 w*tfbb ] 

+ ww* 2 tfb 2 + w 2 w*b(tf) 2 + ww* 2 btfb + ww* 2 b 2 b ] ]i) 

{wtf + w*b) 4 ip = [w 3 w*({tf) 2 btf + 6 t 6(6 t ) 2 ) 

+ w 2 w* 2 ((b f ) 2 b 2 + {tfbf + bW + b{tf) 2 b + (66+) 2 

+ b 2 {tf) 2 ) + w* 3 w(bb*<b 2 + b 2 tfb)]tp 

(wtf + w*bftp = [w 3 w* 2 {(tf)Hb% + {tf) 2 b 2 tf + 

b(tf) 2 b + b\btf) 2 + & t 6 2 (fe t ) 2 + 6(6 t ) 2 66 t + btfb{b ] ) 2 ) 

+ w 2 w* 3 {{tfb) 2 b + 6 + 6 2 6 t 6 + b(tf) 2 b 2 

+ b{tfb) 2 + 6&W + b 2 {tf) 2 b + b{btf) 2 )]ip 



expressed in terms of the creation and annihilation op- 
erators of the bath, b and b^ and do not contain powers 
higher than 3 in their respective expansions and gener- 
ate as many linearly independent vectors as possible on 
T^(Ai), while operating on the state £. With the above 
linearly independent vectors we could construct the new 
control system given by, 



dt 



E + X>* 6 & !«*)> +E ff i J) (.9& t +g*b)m)) +£ Mll a«( W 6t + w * b y\m) 
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i=0 



j2u 2i aW( W tf +w*by\n(t)) +$>3*4 2) M> t +™*by\m) +E^4 2) M> t +^m*)> (47) 

I 



i=0 



i=0 



For the control system described above where the control entire control algebra and hence, 

vector fields {Kjt}, < i < 5 and 1 < j < 4, span the [A) R ^ c A + ^ < i < 5 and 1 < j < 4 (48) 
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where Q = {K lt ■■■ , K 24 } LA = span{ifi, • • • , K 24 }. It 
now remains to know if there exists a controlled invariant 
distribution A, that satisfies the condition stated above. 
It can be seen that since A is a subspace of the tangent 
space T^(Ai) at £, the equation above is trivially satisfied. 
The only additional constraint that A needs to satisfy is 
that it be a part of the ker(dy), the nullspace of y at the 
point £, which is a subspace of the tangent space T^(j\4) 
itself. kei(dy) is comprised of vectors of the form -ff |£) 
where if is a real linear combination(with coefficients 
possibly a function of the state £) of skew hermitian op- 
erators, with the additional constraint that, L^^y = 0, 
which translates to the commutator, [C, H] = 0. Since 
the covector dy is one dimensional for a scalar function y, 
the corresponding nullspace ker(dy), would be n — 1 di- 
mensional where n is the dimension of the tangent space. 
Some of the vectors in ker(dy) are, 

(iW ^ if 2 ))^ + w*by\£(t)), 1 < i < 5 
(of) + aW)(wtf + w*b) l \£(t)), l<i<5 

+ aW)(wtf - w*bY\£(t)), 1 < i < 5 etc 

where the operators a z , 1, above are to understood as 
skew hermitian operators as before. It is to be noted that 
the algorithm presented in the previous section would 
terminate after the first iteration as the condition is al- 
ready satisfied and would yield ker(dy) as A*, the max- 
imal controlled invariant distribution. The least value 
that A could take according to the necessary conditions 
of theorem IVIII.ll Ki G A C ker(dy) is, the one di- 
mensional vector space spanji'Tj}, itself. The algorithm 
presented in the previous section is designed to yield the 
maximal invariant subspace, which guarantees decoupla- 
bility. But in order to compute the feedback we could 
work with any A that is a subspace of maximal A* and 
contains the minimal spanjif/}, as long as the condition 
(gS]| is satisfied. 

Feedback Synthesis In order to determine the feed- 
back let us work with the minimal A = spanji^/}. 
It is possible to construct n — 1 vectors where, n = 
2 x dim(T^(M)) vectors v 2 ,---v n S T^(M) that com- 
mutes with v\ = Ki, (i.e) [i>i,i>j] = 0. Reindexing the 
control vector fields as K\ ■ ■ ■ K r , where r = n = 24 in 
this case, and since Ki span the tangent space we can 
write, 



(49) 



i=l 



where d is a non-singular real matrix. Hence we could 
rewrite, 



fK 1 \ 

K 2 



\K r ) 



dn ■ ■ ■ dir 



drl 



fv 1 \ 



i' 2 



\v r J 



Following the proof of Theorem IVII.ll we can form the 
vectors, K° = S.v, where, S — d^ 1 but with first column 
replaced by zeros, (i,e), 



S = 



S12 • • • Sir 



Now, the feedback parameter (3 is such that 



(3 x 



K° 2 

\k°) 



( V2 \ 

V3 



v r 

V o J 



/ 1 ■ • ■ \ 
1 0--- 



••• 1 

Vo o ••• o o/ 



v 2 



,\v r J 



J 



Since the above equations holds for all v and K vectors 
we could write f3(£,)S = J, but since the above equation 
remains unaltered when S is replaced with d -1 , we can 
calculate the feedback parameter as, (3 = J.d. The closed 
loop vector fields are given by K = (3(^)K . Similarly the 
parameter a can be calculated by incorporating Kq in 
the equation. For any Kq — J2i=i c j v ji we can ^ n< ^ an & 
such that, 



K 



CiVi 



(50) 



for some c\ as a function of the state £. The parameter 
a is given by a — a. [3 and the closed loop drift vector 
field is given by, Kq = ^ oii.Ki. It can be seen that the 
above closed loop vector fields as in the proof does satisfy 
invariance w.r.t A, (i.e) [A, Ki] C A,0 < i < r. Hence 
the system is completely decoupled even in the presence 
of symmetry breaking control hamiltonians via classical 
state feedback. 

**we cannot find a suitable basis transformation using 
real matrices to a known set of commuting vectors such 
as ci|000),ci|102), • • • <E T^M) etc where ci,c 2 € C, as 
performed in where vectors were transformed to co- 
ordinate basis in W 1 in order to determine the feedback. 
Hence the task of finding commuting vectors were simpli- 
fied by such a transformation in the classical case. The 
difficulty is due to fact that (i) coefficients of the states 
complex and effectively carry twice the dimension, (ii) 
tangent vectors at point £ is different from that of an- 
other point £i, hence a fixed coordinate transformation 
does not work for every £. Whereas in the case of R™ 
tangent space at every point x is the same. 

It can also be noted the controllability properties of 
the system are unaltered in the presence of feedback. 
The problem of disturbance decoupling is that of modify- 
ing the observability of the control system via feedback. 
It is very well known from classical control theory that 
feedback can modify the observability properties of any 
system but not the controllability properties. It is the 
observability of the decoherence that we intend to mod- 
ify in the above work by modeling it as a disturbance 
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decoupling problem thus rendering the decoherence act- 
ing on the system unobservable on the states of interest. 
However in order to accomplish the goals we had to in- 
troduce additional couplings and a bait subsystem that 
were not a part of the system initially. 



quantum noise is tensorial. The noise parameter g^ and 
gl dictate the coupling between the environment and the 

system , (i.e), K 1 = (aP + a? ] ) ® (g* k b k + g k b\)\£) cor- 
responds to the classical noise vector p(x), and it can be 
easily seen that there is no noise operating on the system 
in the classical sense. Hence decoherence is not classical 



XI. INTERNAL MODEL PRINCIPLE 

In order to decouple the output from the environment 
one needs to determine the feedback coefficients and 
/3(£) where both depend on the combined state of the sys- 
tem and environment. Hence one needs to have a good 
estimate of the system as well as the environment itself 
for successful implementation of feedback decoupling. In 
other words the state observer must include a model of 
the environment which would enable us estimate its state. 
At this point, the important differences between classical 
and quantum decoupling problems can be understood at 
the outset. The necessary condition in terms of the oper- 
ator algebra [C,Hsb] C C was instrumental in design of 
the bait subsystem. However the structure of the system 
needed to be altered in order to, 

(i) Artificially induce coupling between qubits 1, 2 and 
the environment with the help of the bait. 
(ii) Generate vector fields in higher power of the envi- 
ronment operator to as to generate linearly independent 
vectors. 

Hence it was necessary to modify the core system in 
more ways than one in order to perform decoupling. 
Hence, even though environment is an undesirable in- 
teraction the higher powers of the same helped us gen- 
erate linearly independent vectors in the tangent space, 
which was absolutely necessary for decoupling. Hence 
the environmental coupling here befits the description of 
necessary evil. In classical dynamic feedback Q the de- 
sign of controller depends on the exosystem. In contrast 
the state observer/estimator needs to know the model of 
environment in order to estimate the combined state £ 
and calculate the feedback. Hence the model discussed 
above could be thought of as the Internal Model Princi- 
ple analog of quantum control systems. In addition clas- 
sical output regulation problem concerns with following 
a reference signal in the presence of environmental dis- 
turbace that depends on a prescribed exosystem. On the 
other hand the disturbance decoupling problem focusses 
on eliminating the effects of the environment. 



(ii) Vector spaces in quantum control systems are over 
complex fields. This increases the dimensionality by 2 
fold in many instances where linearly combination has 
to be taken. Hence in order to generate every vector in 
a vector space of n independent states, we require 2n 
linearly independent vectors. 

(iii) The necessary and sufficient conditions impose re- 
strictions on the form of control hamiltonian that could 
help decouple the system. From the conditions derived 
above, it is impossible to decouple one part of the sys- 
tem from the other unless our control hamiltonians op- 
erate on the both the hilbert spaces non-trivially (i.e) 
Hi e B{TLa®'Hb), the set of linear operators in the joint 
hilbert space of both the systems. It was in light of this 
condition that the bait system was originally introduced. 

(iv) Distributions need not necessarily be singular. For 
instance the tangent space of an su(2) system is spanned 
by or z \^),a x \£_),o- y \0^\0, where £ = c |0) + a\l) and 
the operators are again assumed to be skew hermitian 
counterparts of hermitian o~ z ,o~ x ,o~ y . Even though the 
four vectors are linearly independent for almost all non- 
zero values of cq and C\ the distribution is non-singular. 
Consider |£) = |0) and the corresponding tangent vectors 
are — i|0), i|l), — i|0), whose real linear combination is 
rank deficient. Hence it can be seen that the vector |0) 
does not belong to tangent space T| ) at the point |£) = 
|0). In general the tangent vectors at point £ is different 
from that of another point £i. One of the most serious 
implications is that we cannot find a linear map that 
transforms the distribution A to a constant d dimensional 
distribution, 



T.A = 



jdxd 





at every point £, an approach that was used in 
Isidori[J| to greatly simplify finding commuting vectors 
|t>i), • • • \v n ) in an n dimensional tangent space. The com- 
muting vectors were just taken to be the co-ordinate basis 
at every point x. 



XII. BILINEAR INPUT AFFINE 
REPRESENTATION OF QUANTUM SYSTEMS 

In this section we will attempt to highlight a few more 
important differences between the decoherence control in 
quantum systems and disturbance decoupling of classical 
input affine systems in R n . 

(i) Classical noise is additive, x — f(x) + Uig^x) + 
wp(x) and operate on the same vector space. Whereas 



XIII. CONCLUSION 

In this work we provided the conditions and a step 
by step procedure to calculate a classical deterministic 
feedback under which the 2-qubit system could be suc- 
cessfully decoupled from decoherence. As mentioned be- 
fore the analysis carried out in the bilinear form only 
helped us learn about the control hamiltonians helpful in 
decoupling the system but also provided a solution un- 
der which the system would be completely decoupled as 
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opposed to partial or nth order decoupling discussed in 
various previous work. Such a control strategy would be 
immensely helpful in performing decoherence free quan- 
tum computation thus enabling us to exploit the com- 



putational speed up provided by quantum parallelism. 
However in order to determine the feedback one needs to 
have a good estimate of the state of the system. 
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